Nonequilibrium spin fluctuations in single-electron transistors 
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We show that nonequihbrium spin fluctuations significantly influence the electronic transport 
in a single-electron transistor, when the spin relaxation on the island is slow compared to other 
relaxation processes, and when size effects play a role. To describe spin fluctuations we generalize 
the 'orthodox' tunneling theory to take into account the electron spin, and show that the transition 
between consecutive charge states can occur via a high-spin state. This significantly modifies the 
shape of Coulomb steps and gives rise to additional resonances at low temperatures. Recently 
some of our predictions were confirmed by Fujisawa et al. [Phys. Rev. Lett. 88, 236802 (2002)], 
who demonstrated experimentally the importance of nonequilibrum spin fluctuations in transport 
through quantum dots. 

PACS numbers: 73.23.Hk, 72.25.-b, 73.40. Gk 
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Electronic transport in single electron transistors 
(SETs) with metallic islands or semiconducting quan- 
tum dots (QDs) was usually described in terms of the 
'orthodox' theory of sequential tunneling [1]. The relax- 
ation processes on the island were characterized only by 
the energy relaxation time t^, with the tacit assumption 
that the spin-flip relaxation time Tgf is short or compa- 
rable to Te [2-4]. In real systems, however, the ratio 
Tgf/Te can be as high as 10^ - 10'' [5-7]. In this letter 
we show that, when the spin relaxation time is signif- 
icantly longer than the energy relaxation time Te and 
injection time ti, large nonequilibrium spin fluctuations 
(NSF) arise on the island. When the density of states 
for the island is low, e.g., due to size effects, the spin 
fluctuations lead to fluctuations of the spin-splitting of 
the chemical potential. This, in turn, can significantly 
modify the transport characteristics. In particular, the 
Coulomb steps are smeared because transitions to higher 
charge states may occur via a set of associated spin states, 
each of them corresponding to different temporary elec- 
trochemical potentials. Such NSF can limit the perfor- 
mance of spin electronic devices. The importance of spin 
fluctuations has been recently conflrmed experimentally 
by Fujisawa et al. [8] for a few-electron QD. 

The effects described here differ from other phenom- 
ena related to energy level quantization and electron spin, 
which have been investigated recently theoretically [2,3] 
and experimentally for paramagnetic [9-11] and ferro- 
magnetic grains [12,13]. They also differ from further 
spin-related phenomena, like spin blockade due to ex- 
change interaction [14], parity effects [15], Kondo effects 
[16,17], and quantum entanglement. 

Before we analyze spin effects in normal metal SETs, 
we consider first the corresponding double tunnel junc- 
tion without Coulomb blockade and with a linear current- 



voltage relation. Tunneling of a single electron increases 
or decreases the magnetization M of the island by 1. 
Here M 



Nf 



Ni, where iV| 



(Ni) are the num- 
bers of excess electrons with spin a =t (i). In the 
absence of intrinsic spin-flip processes, the time evolu- 
tion of the magnetization M{t) can be mapped onto 
a one-dimensional diffusion process with (M) = and 
(A/2) 1/2 ^ There are, however, two main processes 
which restrict the increase of the fluctuations: (i) spin- 
flip relaxation processes, and (ii) spin splitting of the 
electrochemical potential /i| ^ fi^ (two spin components 
of the island can be treated as two independent elec- 
tron reservoirs with different electrochemical potentials 
fia) due to spin accumulation M [18], which modifies the 
tunneling probability. In the non-equilibrium situation 
due to an applied bias V, fluctuating numbers N^- give 
rise to nonequilibrium fluctuations Na-AE of the electro- 
chemical potentials Ha- (here AE = 1/D{Ey) and D{Ef) 
denotes the density of states at the Fermi level in the 
island), which influence the transport current. The cur- 
rents flowing through the left L junction for both spin 
orientations are /lt = 1/2 (VRi^/R -\- MAE/2e)/Ri, 
and hi = 1/2 {VRi^/R - MAE/2e)/Ri^. i?R (i?L) 
denote the resistance of the right (left) junctions and 
R = i?R -I- i?L- However, the total current /l = /lt + 
/lx = V/R is independent of M. Straightforward cal- 
culation shows that for eV/AE ^ 1 and for low tem- 
peratures the magnetization fluctuations are described 
by a Gaussian distribution with the standard deviation 
(M2)1/2 « y'ieV/AE) a/(a + 1)2, where a = Ru/Rl- 

In contrast to the example discussed above, NSF in sys- 
tems with nonlinear current-voltage characteristics can 
influence the average transport current. We will describe 
now how spin fluctuations influence transport through 
a SET, where the charge fluctuations are strongly sup- 
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pressed due to the Coulomb interaction, but spin fluc- 
tuations can be relatively large. For charging effects 
the relevant charging energy is Eq = e^/2C, where 
C = Cl + Cr + Cg is the total capacitance of the island, 
which is the sum of capacitances of the left (Cl) and 
right (Cr) junctions and of the gate (Cg). The electronic 
transport through the system in a stationary state is gov- 
erned by the solution of the generalized master equation 
[19] 

0=-{T{N^,Ni)+nu{N^,N^) + ni,^{N^,Ni)}PiN^,Ni) 
+T+{N^-1,N^)P{N^ -1,N^)+T+{N^,N^ -1)P{N^,N^ -1) 
+rf (AT^ +l,Ni)P{N^ +l,Ni)+ri{N^ +l)P{N^,Ni +1) 
+n^,iiN^ - 1, iVi + l)-P(iVT - 1, iVi + 1) 

+ 1, - i)^W + 1, - 1) (1) 

for the probability P{N'^,Ni) to find iV| and Ni excess 
electrons on the island ( A'^ = N-^+N^ is the total number 
of excess electrons). The first term in Eq. (1) describes 
the rate at which the probability of a given configura- 
tion decays due to electron tunneling to or from the is- 
land, whereas other terms describe the rate at which this 
probability increases. The terms account for spin-fiip 
relaxation processes. The coefficients entering Eq. (1) 
are defined as r±(7V|,7V^) = E.=LRr^<.(A^T.^i) and 
m^Ni) = Ep=±E.rP(7VT,iVi), where r± (7VT,iVx) 
are the tunneling rates for electrons with spin a, tunnel- 
ing to (+) the grain from the lead r = L, R or back (— ). 
These coefficients are given by the following expression: 

Tt,{N^,N^)=J2lLF^{Eia\N^,N^)f^{Ei,+E+iN)-EF), 

i 

i 3 

X F^{E^^\N^,N^). (2) 

Here, f^{E) (/~ = 1 - /+) is the Fermi distribution 
fimction, whereas F+(£,^|iV|, ]V|) [20] (F^ = 1 - F+) 
describes the probability that the energy level Ei^ is oc- 
cupied by an electron with spin a for the particular con- 
figuration (iV-f,7V|). The parameter is the tunneling 
rate of electrons between the lead r and the energy level 
Ei^ of the island, and coiaja is the transition probability 
from the state ia to ja due to the spin-fiip processes. 
The energies (^) and E^{N) are given by E^{N) = 
Cn/C eV + C/±(iV) and E^iN) = -Cl/C eV + C/±(iV) 
where C/±(Af) = Ec[2{N - N^) ± 1] and A^^ = C^VJe, 
with Vg denoting the gate voltage. 

From the solution P{Ni,Ni) of the master equation 
(1) we obtain the current flowing through the island 

^- = «E E PiN^,Ni){TUN^,N^)-T;,{N^,N^)} . (3) 

a- Nf,Ni 

In our calculations we assume that discrete energy lev- 
els Eicr are equally separated, with the level spacing AE. 



The tunneling rates arc then given by the formula 
jI^ = AE/e'^Rra [3], where Rra is the resistance of the 
junction r for spin a {Rri = Rri = 2Rr, in the case 
under consideration). For spin- flip processes we assume 

In Fig. 1 we show calculated transport characteris- 
tics for symmetric (left column) and asymmetric (right 
column) junctions. In both cases the transport char- 
acteristics are calculated in the fast (ti ^ Tsf,Te) and 
slow (Tst 3> Ti ^ Te) spin-fiip relaxation limits, where 
n ~ l/7f<7 ~ e// (for QDs, it is shown in Ref. [6] that 
Tsf > 1/tis, Te < Ins and rr ~ 1 10ns for 7 ~ 0.1 InA). 
In Fig. l(a,f) the conductance- voltage characteristics are 
presented for both limits. The effect on the current is 
relatively weak (a few per-cent at maximum), as shown 
in Fig. l(b.g). However, the effect on the differential 
conductance can be significantly larger, about a few tens 
per-cent at maximum, as shown in Fig. l(c,h). Generally, 
the most pronomiced effect of the spin fluctuations is on 
the nonlinear parts of the transport characteristics. This 
effect is similar to the one produced by an increase in 
effective temperature Tcs [22,23] of the system (different 
from the bath temperature T). However, both effects 
can be easily distinguished because of the parabolic de- 
pendence on the bias voltage V of the former effect (see 
the discussion below). For intermediate spin-flip times, 
0.1 < Tsf/ri < 10, the conductance smoothly crosses over 
from one limit to the other. 

In Fig. l(d,i) we show the bias dependence of the spin 
fluctuations (M^)^/^ in the limits of both short and long 
spin relaxation times. In the former case the fluctua- 
tions are almost constant and close to 1. In the lat- 
ter case and for symmetrical junctions, the fluctuations 
nearly follow the law for junctions without Coulomb ef- 
fects ((M2)1/2 ^ ^{eV/AE) a/(a+ 1)2). However, for 
asymmetric junctions there are pronounced oscillations in 
the fluctuations amplitude with increasing bias V , which 
are related to the charge accumulation and correlated 
with the charge fluctuations [21] (see Fig. l(j)). For com- 
parison, we show in Fig. l(e,j) the amplitude of charge 
fluctuations. For asymmetric junctions the oscillations 
are due to charge accumulation. 

In the inset of Fig. 1(e), we show the probability P of 

the = 1 state in the transition range (from A'' = to 
A'^ = 1) for the symmetrical junction. One can see that 
NSF in the case of slow spin-flip relaxation allow the 
transition to the next charge state at lower bias voltage. 

Figure 2 demonstrates the mechanism how the spin 
fluctuations assist the system to enter the higher charge 
states. For strong NSF, an electron needs a lower energy 
to enter the island, because some of the double degen- 
erate states below the Fermi energy are empty (com- 
pare Fig. 2(b) and Fig. 2(a)). Figure 2(c) illustrates 
that the onset of the transition to the next charge state 
is linked to the existence of a high spin state, and the 
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new charge state at the onset can appear only via the 
high spin state. The effect is relevant as long as the 
energy of thermal fluctuations is lower than the energy 
of NSF, 2(Af2)i/2A£; > k^T. Thus, even if the dis- 
crete states are not resolved {k^T > AE'), the effects 
due to spin fluctuations can be important. In Fig. 3 
we show the amplitude of NSF as a function of Ai?. 
From this figure follows that the NSF effects can be 
observed in electronic transport when AE/Eq ^ 0.01. 
This condition can be easily achieved for typical QDs 
[1] of radius 200nm, 2Ec ~ ImeV, AE ~ O.OSmeV and 
for T < IK. For smaller QDs the effect is important 
also at higher temperatures. For larger islands with high 
density of states D{E-p) the fluctuations of M have very 
large amplitude -^/eV/A^), but the associated effec- 
tive spin splitting of the electrochemical potential MAE 
is small (~ V eV AE). For metallic systems the condition 
2(M2)V2a£: > knT could be fulfllled for small metal- 
lic particles with diameter in the range of several nm. 
For strong exchange interactions QDs can have a ground 
state with spin S > 1/2 [14,24]. But even in this case 
spin fluctuations are important, as pointed out by Kleff 
and von Delft [13] for ferromagnetic grains. 

The time Tgf depends mainly on the strength of the 
spin-orbit interaction and on the impurity contents [7]. 
Thus counter-intuitively, for SETs with dirty island and 
strong spin scattering one should expect sharper trans- 
port characteristics. The NSF effects can be detected by 
measurements of the conductance step or peak widths as 
a function of temperature. The widths should saturate 
at a value that scales with the bias voltage as V^/"^. In 
the absence of NSF there will be no saturation observed, 
or the saturation will be at a value (related to some T^s) 
which does not depend on the bias as y^/^. 

Very recently, Fujisawa et al. [8] obtained results for 
very small QDs at low temperatures, which confirm the 
importance of NSF and which can he accounted for by 
our model. In Fig. 2(d) we show schematically the po- 
sitions of the conductance peaks (current steps), i.e., dc 
excitations in the low-temperature limit, k^T ^ AE, 
where discrete energy levels are resolved. In addition to 
the standard resonance peaks (corresponding to the dot- 
ted lines), there are new ones indicated by the dashed 
lines which start inside the diamond corresponding to 
single electron tunneling transport (SETT) (not at its 
border) when a particular spin excitation appears. Dou- 
ble electron tunneling transport (DETT) can occur al- 
ready within the SETT diamond region due to NSF (the 
regions marked in gray in Fig. 2(d), where also the mag- 
netic state of the excitations is indicated). We find that 
not only in results of Ref. [8] but also in earlier data (see 
e.g. Fig.l.(A) in Ref. [25]) the borders between SETT 
and DETT diamonds could be affected by spin fluctua- 
tion. The new resonance peaks are related to the fact 
that the ratio Eq/AE is generally not integer, and the 
electron levels for adjacent charge states are effectively 



shifted by mod{Ec,AE). 

One may expect a similar effect of NSF on the cotun- 
neling current, too. The fluctuations can influence in 
particular the first step in the current- voltage character- 
istic (the first peak in the conductance). One may also 
expect that NSF modify the current shot noise. 

In summary we have generalized earlier spinless de- 
scriptions of SETs by taking into account the fact that 
spin-flip relaxation time Tgf is usually much longer than 
the energy relaxation time r^. The spinless models de- 
scribe properly only the situation when Tgf is short or 
comparable to r^. When, however, Tgf » r^, large 
nonequilibrium spin fluctuations may occur in the is- 
land. These fluctuations have a significant influence on 
the transport characteristics. Our model accounts for 
recent experimental observations [8]. 
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FIG. 1. (a,f) Differential conductance vs. bias voltage 
V for two limits: short (dashed line) and long (solid line) 

spin-flip relaxation time; (b,g) relative difference (/i — /s)/7s, 
whore Ji and Is arc the currents for the two limits; (c,h) rela- 
tive difference (Gi — Gs)/Gs; (d,i) spin fluctuations (M^)^/^; 
(e,j) charge fluctuations ((A^^) - {Nff^. The curves are 
calculated for symmetric (left column) and asymmetric junc- 
tions (right column). Inset: The probability P of the 
N = 1 charge state for a symmetrical junction as a function 
of the bias voltage V. For symmetric junctions Cr = Cl, 
Rr = i?L and the gate voltage = 0.5 + 1/4 AE/Ec (at 
resonance), whereas for asymmetrical junctions Cn/Ci, = 3, 
Rr/Rl = 100 and the gate voltage A^^, = 0. The other pa- 
rameters are ksT = OAAE and AE/Ec = 0.1 for both cases. 
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FIG. 2. Examples of specific spin configurations of the is- 
land for the same charge state of the system with (a) short 
and (b) long spin-flip relaxation times. Part (c) shows how 
the system enters the next charge state. The probability 
P(iV|,iV|) for N = (iVi+A^i) = 1 is shown there as a function 
of M = (A''! —Ni) for several values of the bias voltage V. All 
parameters are the same as in Fig. 1(a). The lines are a guide 
to eye only. Data for different V arc offset vertically. (d) 
The scheme of conductance peaks (current / plateaus) for dc 
excitation transport in the V — Vg plane for symmetric junc- 
tion, Ec/AE = 4.5 and ksT < AE. Sohd lines determine 
the Coulomb blockade (CB) diamond structure indicating also 
the onset of single electron tunneling transport (SETT) and 
double electron tunneling transport (DETT) (two excess elec- 
trons on the island are possible) in the absence of NSF. Dotted 
lines are usual efiiects related to offset in transport of the next 
discrete energy levels. Dashed and dotted-dashed lines indi- 
cate the onset of the consecutive charge states due to particu- 
lar spin excitations (as indicated) which without spin fluctua- 
tions are not accessible. Dashed lines indicate also formation 
of new resonances due to NSF. Symmetry is broken due to 
the parity effect. 
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FIG. 3. Spin fluctuations (M^)^/^ (dashed line) and 

the corresponding splitting of electrochemical potential 
{mY^'^AE/Ec (solid line) as a function of energy level spac- 
ing AE/Ec, calculated for ksT = OMEc and eV/2Ec = 3. 
The other parameters are as in Fig. 1(a). 
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